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Abstract 

The aim of this paper is to prove central limit theorems with respect to the annealed measure for 
the magnetization rescaled by '/N of Ising models on random graphs. More precisely, we consider the 
general rank-1 inhomogeneous random graph (or generalized random graph), the 2-regular configuration 
model and the configuration model with degrees 1 and 2. For the generalized random graph, we first 
show the existence of a finite annealed inverse critical temperature 0 < /3c" < oo and then prove our 
results in the uniqueness regime, i.e., the values of inverse temperature j3 and external magnetic field 
B for which either /3 < /3“ and R = 0, or /3 > 0 and B ^ 0. 

In the case of the configuration model, the central limit theorem holds in the whole region of the 
parameters P and B, because phase transitions do not exist for these systems as they are closely related 
to one-dimensional Ising models. Our proofs are based on explicit computations that are possible since 
the Ising model on the generalized random graph in the annealed setting is reduced to an inhomogeneous 
Curie-Weiss model, while the analysis of the configuration model with degrees only taking values 1 and 
2 relies on that of the classical one-dimensional Ising model. 


1 Introduction and main results 

1.1 Motivation 

1.1.1 Ising models on random graphs 

The ferromagnetic Ising model is the most well-known example of statistical mechanics system describing 
cooperative behavior. Its probabilistic formulation |18) amounts to an infinite family of random variables 
taking values in { — 1,1} (so-called spins) whose joint law is given by the Boltzmann-Gibbs distribution. 
The properties of such families of random variables are crucially determined by the spatial structure where 
the spin variables are sitting. For instance, for the Ising model on with nearest-neighbor interactions, 
the model displays a second-order phase transition for d > 2. Furthermore, the universality prediction 
states that the precise details of the interactions are not relevant for the near-critical behavior, so that 
around the critical temperature each universality class is described by a single set of critical exponents. 

Besides regular lattices, in recent years much attention has been devoted to the setting in which the 
spin variables are placed on the vertices of random graphs liiiiniiiiiiiaiisiiiiiiisiiiiiiiTiEaiEsj. Such 
random graphs aim to model emergent properties of complex systems consisting of many interacting agents 
described by a network. Several studies on empirical networks have found that two random elements of 
the network are typically within relatively short graph distance (the so-called small-world paradigm). 
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whereas there does not exist a typical scale for the number of neighbors that a random element has (the 
so-called scale-free paradigm where degrees in the network are proposed to have a power-law distribution) 
[261 [271 da [2D]. As a consequence, there are vertices with very high degree that often play an important 
role in the functionality of the network. 

Thus, the combination of the ferromagnetic Ising model on a random graph describes situations in 
which single units establish macroscopic cooperative behavior in the presence of the random and complex 
connectivity structure described by a network. Here the playing field has two levels of randomness: hrstly, 
the probabilistic law of the spins and, secondly, the probability distribution of the graph. So far, most of 
the studies have focused on the so-called (random) quenched state, in which the random graph is considered 
to be hxed once and for all. In this paper, we instead consider the annealed state: the Ising model at every 
time sees an average of the possible random graphs [D] [22], rather than one realization of the graph. The 
annealed measure is particularly relevant for applications in socio-economic systems, in which the graph 
dynamic models the evolution of social acquaintances, or the brain, in which graph edge rearrangements 
represent the evolution of synaptic connections. We explain the role of the annealed and the quenched 
laws in more detail in the following section. 

1.1.2 Annealing 

To understand the difference between the quenched and annealed settings, it is convenient to think of 
a microscopic dynamics yielding the equilibrium state. For instance, one could imagine that the spins 
are subject to a Glauber dynamics with a reversible Boltzmann-Gibbs distribution and the graph also 
has its own dynamical evolution approaching the graph’s stationary distribution. In general, these two 
dynamics are intertwined and both concur to determine the equilibrium state, i.e., the asymptotic value 
of an ergodic dynamical time average. The quenched and annealed state arise as follows: 

(a) In the quenched state, the changes of the graph happen on a time-scale that is infinitely longer 
than the time-scale over which the changes of the spin variables occur. Thus in the quenched state 
the graph viewed by the evolving spins is frozen. One distinguishes between the random quenched 
measure, i.e., the random Boltzmann-Gibbs distribution of a given realization of the graph, and the 
averaged quenched measure, i.e., the average of the Boltzmann-Gibbs distribution over the graph 
ensemble. Several thermodynamic observables (e.g., the free energy per particle, the internal energy 
per particle, etc.) are self-averaging, and therefore the random quenched values and their averaged 
quenched expectations do coincide in the thermodynamic limit. In the study of the fluctuations of 
the properly rescaled magnetization one finds a Gaussian limiting law. Interestingly, the asymptotic 
variances of the random quenched and averaged quenched state might be different [8] due to local 
Gaussian fluctuations of graph properties. 

(b) In the annealed state, the environment seen by the spins includes all possible arrangements of the 
random graph. The annealed measure (defined later in (11.71) 1 is given by the stationary reversible 
measure of a Glauber spin dynamics in which the transition from a conhguration cr to another 
conhguration a' occurs with probability 


]E[e-/3//(a')] 


( 1 . 1 ) 


where H is the Hamiltonian and E[-] represents the average over the graph ensemble. The above 
dynamics corresponds to an extremely fast random graph dynamics in which we do not even observe 
the graph at any time, but merely see it averaged over the random graph distribution. This is 
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equivalent to an effective Glauber dynamics with (annealed) Hamiltonian equal to 

//“(a) = log(E[e-^^(‘")]). (1.2) 

Thus, by construction, the annealed measure is necessarily non-random. We will be interested in 
the properties of the Gibbs measure corresponding to the dynamics in (II.ip . which corresponds to 
the stationary or infinite-time distribution of the spins under the dynamics. While the Glauber 
dynamics m corresponds to infinitely fast graph dynamics compared to the spin dynamics, the 
stationary distribution can equally well be viewed as a dynamics where the graph and the spin evolve 
at equal speeds, as is the more usual viewpoint in statistical mechanics. Note that, in the definition 
of the annealed pressure (see (jl.Sp ). the averages taken w.r.t. the spins and the graph are completely 
symmetric, which can be seen as another argument in favor of the view that the corresponding 
dynamics run equally fast and that the limiting measure corresponds to the average w.r.t. graph 
and spins alike. In this paper, we will study annealed central limit theorems for the ferromagnetic 
Ising model on random graphs, in order to deduce what the effect of annealing on the macroscopic 
properties of the Ising model is. 

The definition of the annealed measure in the context of Ising models on random graphs is thus different 
than in other class of problems with disorder, such as random walks in random environment [5]. In that 
context, annealing is rather similar to what here we have called the averaged quenched measure. 

In disordered systems (such as spin glasses |241 [6]). annealed disorder is usually considered to be easier 
to deal with mathematically, since the average on the disorder and the thermal average are treated on 
the same footing. This is true whenever the edges of the graph are independent, due to the form of the 
Hamiltonian that allows a factorization of expectations w.r.t. the bond variables. If instead the edge 
distribution in the graph does not have a product structure, the annealed case can actually be more 
difficult than the quenched case. Indeed, whereas the random-quenched case is dominated by the typical 
realization of the graph (often having the local structure of a random tree), in the annealed case (as in 
the averaged-quenched case) the rare graph samples actually give a contribution that can not be ignored. 
This is due to the fact that the Ising model gives rise to exponential functionals on the random graph, 
and expectations of exponential functionals tend to be dominated by rare events in which the exponential 
functional is larger than it would be under the quenched law. Deriving such statement rigorously requires 
a deep understanding of the large deviation properties of random graphs, a highly interesting but also 
challenging topic. 

In this paper, we consider graph ensembles of both types, i.e., random graphs with independent edges 
(these are generalized random graphs) or dependent edges (in this case, we study the configuration model). 
These are described in the following section. 


1.2 Random graph models 

We denote by a random graph with vertex set = [N] and edge set C x V^. 

Here and in the rest of this paper, we write [N] = {I,..., N} for the vertex set of G^- For any iV E N, 
we denote by Qjv the probability law of the random graph Gff. In this work, we consider two classes 
of random graphs: the configuration model and the generalized random graph. We next introduce these 
models. 
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1.2.1 The Generalized Random Graph 


In the generalized random graph, each vertex i G [A^] receives a weight Wi > 0. Given the weights, edges 
are present independently, but the occupation probabilities for different edges are not identical, instead, 
they are moderated by the vertex weights. For a given sequence of weights w = the graph is 

denoted by GRGjv('if)- We call the Bernoulli indicator that the edge between vertex i and vertex j is 
present and pij = P {lij = 1) is equal to 


WiWj 

Pij = -n —^^, 
tjv + WiWj 

where t'jv is the total vertex weight given by 

N 

£n = '^Wi . 

i=l 


(1.3) 


(1.4) 


Denote by the weight of a uniformly chosen vertex Ijv G [A^]. The weight sequence of the 

generalized random graph GRGjv('U^) is often assumed to satisfy a regularity condition, which is expressed 
as follows: 


Condition 1.1 (Weight regularity). There exists a random variable W such that, as N ^ oo, 

(a) A W, 

(b) E[VFjv] = ^ Eie[JV] ^ < OO' 

(c) E[W2] = ^ Eie[iv] < oo. 

T> 

where —> denotes convergence in distribution. Further, we assume that E[VF] > 0. 

In the following, we will consider deterministic sequences of weights that satisfy Gondition 11.11 In 
many cases, one could also work with weights w = (rCi)ig[Ar] that are i.i.d. random variables. For the 
annealed setting, however, one has to be careful, as we will argue in more detail in Section 11.5.11 below. 
Indeed, when the weights are themselves random variables, they introduce a double randomness in the 
random graphs: firstly there is the randomness introduced by the weights, and secondly there is the 
randomness introduced by the edge occupation statuses, which are conditionally independent given the 
weights. Whereas the thermodynamic properties (pressure, magnetization, etc.) in the quenched measures 
of the Ising model on GRGjv('U^) are not affected by the choice of deterministic or random weight sequences, 
the pressure of the annealed Ising model becomes infinite when the weights have sufficiently heavy tails. 


1.2.2 The Configuration Model 

The configuration model is a multigraph, that is, a graph possibly having self-loops and multiple edges 
between pairs of vertices. Fix an integer N and consider a sequence of integers d = (di)ig[Ar]. The aim 
is to construct an undirected multigraph with N vertices, where vertex j has degree dj. We assume that 
dj > 1 for all j G [N] and we denote the total degree in the graph by 

£n •= ^ ^ di- (1-5) 

ie[N] 

We assume to be even in order to be able to construct the graph. 
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Assuming that initially dj half-edges are attached to each vertex j E [A^], one way of obtaining a 
uniform multigraph with the given degree sequence is to pair the half-edges belonging to the different 
vertices in a uniform way. Two half-edges together form an edge, thus creating the edges in the graph. 
To construct the multigraph with degree sequence d, the half-edges are numbered in an arbitrary order 
from 1 to iff. Then we start by randomly connecting the first half-edge with one of the iff — 1 remaining 
half-edges. Once paired, two half-edges form a single edge of the multigraph. We continue the procedure 
of randomly choosing and pairing the half-edges until all half-edges are connected, and call the resulting 
graph the configuration model with degree sequence d, abbreviated as CMjv(d). 

We will consider, in particular, the following models: 

(1) The 2-regular random graph, i.e., the configuration model with di = 2 for all i E [A^], which we 
denote by CMjv(2). 

(2) The configuration model with di E {1,2} for all i E [A^], which we denote by CMjv(l,2). In 
CM^(1, 2), for a given p E [0,1], we have N — [pN\ vertices of degree 1 and [pN\ vertices of degree 

2 . 


1.2.3 Properties of GKGjf{w) and CM^(d) 

The existence of a phase transition in the structural properties of the graph depends on the asymptotic 
degree D, i.e. the weak limit, provided it exists, of the sequence {D]\f)]sf>i where Djv is the degree a 
uniformly chosen vertex J^v £ [A^] in the graph. In order to state this result, we introduce some notation 
that we will frequently rely upon. Let the integer-valued random variable D have distribution P = {pk)k>i, 
i.e., ¥{D = k) = pk, for k > 1. We define the size-biased law p = {pk)k>Q of D by 


Pk = 


{k + l)pk+i 
E[D] 


where the expected value of D is supposed to be finite, and introduce the average value of p by 

E[D{D - 1)] 


V := 


:= ^ kpk 


k>0 


K[D] 


( 1 . 6 ) 


For CM]v(2), the asymptotic degree distribution equals P(Z1 = 2) = 1, while for CMjv(l, 2), the asymptotic 
degree distribution equals P(Z1 = 2) = p,¥{D = 1) = 1 — p. For GRGjv(w) with asymptotic weight 
distribution W, the asymptotic degree is a mixed Poisson random variable Poi(lT) where W appears 
in Condition 11.11 i.e., 


P(T) = k) = E 


.,-w 


If 


see e.g., m Chapter 6]. 

It is well known El EU that the above random graphs have a phase transition in their maximal 
component. Indeed, when > 1 a giant component exists, while for iz < 1 the maximal component has 
o{N) vertices. Here, since the degree distribution for GIiGff{w) is D = Poi(VF), we have that v = , 

because E[T)] = E[W] and E[T)(H—1)] = E[1T^]. Thus, depending on W, a giant component for GRGjv(’*a’) 
may exist, while it does not exist for CM^(2) and CMjv(l,2) since, in these cases, n < 1. In fact, for 
CMjv(2), the connectivity structure is quite interesting and explained in more detail in |21| . 
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1.3 Annealed measure and thermodynamic quantities 

We continue by introducing the ferromagnetic Ising model and the annealed measure. We define them 
on finite graphs with N vertices and then study asymptotic results in the limit N ^ oo. We denote a 
configuration of N spins by a, where a is defined on the vertices of the random graph Gpf whose law is 
Qn- 

In our previous work [8], we have considered two Ising models. The random-quenched measure ^Gn{^) 
coincides with the random Boltzmann-Gibbs distribution, where the randomness is given by the graph 
Giv-The averaged-quenched measure is obtained by averaging the random Boltzmann-Gibbs distri¬ 

bution over all possible random graphs, i.e., Pjv(c’') = QAf(l^Gjv 

In defining the annealed measure, the numerator and denominator of the Boltzmann-Gibbs distribution 
g^GN averaged separately with respect to Q^f, as formalized in the following definition: 

Definition 1.1 (Annealed measure). For spin variables a = (ui, ...,(Tjv) taking values on the space of spin 
configurations fljv = { — 1,1}^, we define the annealed measure by 


PN{a) = 


Qn exp 


PE 




a. 


iUj + BY, 


ielAf] 


(Ji 


Q^{Z^{(5,B)) 


(1.7) 


where 


Zn{P,B) = ^ exp 
crGHiv 


P 


+ B (Tj 
ie[Af] 


is the partition function. Here fi > 0 is the inverse temperature and B ^ M. is the uniform external 
magnetic field. 


In this paper, with a slight abuse of notation, we use the same symbols to denote both a measure and 
the corresponding expectation. Moreover, we remark that the measure defined above depends sensitively 
on the two parameters {f5,B). However, for the sake of notation, we will drop the dependence of the 
measure on these parameters. Sometimes we will use Var^(A) to denote the variance of a random variable 
X with law jjL. 

We now define the thermodynamic quantities with respect to the annealed measure: 

Definition 1.2 (Thermodynamic quantities |101 I14|). For a given N £ 'H, we introduce the following 
thermodynamics quantities in finite volume: 

(i) The annealed pressure is given by 

B) = ^ log {Q^ (Z^ i(3, B))). (1.8) 

(a) The annealed magnetization is given by 

Mr,{P,B) = 


where the total spin is defined as 


S^= y a,. 







(in) The annealed susceptibility equals 

X.{I3,B} := = Varp^ , 

We are interested in the thermodynamic limit of these quantities, i.e., their limits as —>■ oo. In this 
limit, critical phenomena may appear. If A4(/3, B) := limTv^oo B), where B) is the average 

of Spi/N with respect to /rG^(-), Pjv(-) or P;v(') and provided this limit exists, criticality manifests itself in 
the behavior of the spontaneous magnetization defined as Ad(/3,0^) = lim5^o+ -^iP^B). In more detail, 
the critical inverse temperature is defined as 

/3c := inf{/3 > 0: M(/3,0+) > 0}. (1.9) 

and thus, depending on the setting, we can obtain the quenched and annealed critical points denoted by 
and respectively. When 0 < /3c < oo, we say that the system undergoes a phase transition at 
(3 = Pc 

From |15| . we recall that, in the general setting of tree-like random graphs to which our models belong, 
the quenched critical inverse temperature is given by 

/3)1'^ = atanh(l/j/), (1-10) 

where n is defined in (m. Let us remark that, in the quenched setting, since < 1 for both CMjv(2) 
and CMjv(l, 2), from (ll.lOli it follows immediately that /3c“ = oo, which means that there is no quenched 
phase transition in these models. In the annealed setting, we will prove the absence of phase transition for 
CMjv(2) and CMjv(l,2) below. On the contrary, we will see that a critical inverse temperature appears 
for GRGjv('if)- 

1.4 Results 

We focus first on the study of the generalized random graph under the annealed measure, obtaining the 
Strong Law of Large Numbers (SLLN) and the Gentral Limit Theorem (CLT) for the total spin Sff. Then 
we present the results in the annealed setting for the configuration models CM^(2) and GMjv(l,2). 

1.4.1 Results for GRGjv(^e) 

The proofs of the SLLN and GLT for GRGjv('it’) require to investigate the uniqueness regime for GRGjv('ii’)- 
For this, we first investigate the existence of the thermodynamic quantities in the infinite volume limit 
with respect to the annealed law. These results will be obtained in the next theorem. They show, in 
particular, that annealing changes the critical inverse temperature. Indeed, the annealed critical inverse 
temperature /3^“ is strictly smaller than the quenched critical inverse temperature /3c when the latter 
exists. In the statement of the theorem below, we will use the notation for the annealed uniqueness 
regime, i.e., 

:= {(/3, R) : /3 > 0, R / 0 or 0 < /3 < Pf^, R = 0} . 

Theorem 1.1 (Thermodynamic limits for the annealed GRGjv('U))). Let (G'jv) 7 v>i ® sequence of 
GRGjv('ii’) satisfying Condition \l.l\ Then the following conelusions hold: 

(i) For all 0 < P < oo and for all B G M, the annealed pressure exists in the thermodynamic limit 
N ^ oo and is given by 

HP,B) := lim ^jv(/3,R), (1.11) 

N^oo 

its value is given in 
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(a) For all {I3,B) E the magnetization per vertex exists in the limit N —oo, i.e. 

M{P,B) := lim Mj,{P,B). 

N—>-oo 

For i? 7 ^ 0 the limit value M{I3,B) equals M{I3,B) = B) and is given by 


M{/3,B) = E 


tanh 


' sinh (/3) 

E[VF] 


Wz* + B 


where z* = z*{l3,B) is the solution of the fixed-point equation 

2 ; = E 


tanh 


/sinh (/3) 


E[H^] 


Wz + B 


sinh (/3) 

E[Ty] 


w 


and W is the limiting random variable defined in Condition \l.l[ 
(Hi) The spontaneous magnetization is given by 

0 if 


M(fi):= lim M((3,B) = 

B^0+ 


fio if fiiw 


( 1 . 12 ) 


and the annealed critical inverse temperature is 

= asinh (1/zz) , 

where v, defined in (ED, is given by u = E[ll^^]/E[hF] and W is the limiting random variable 
introduced in Condition \l.l\ In particular, if v > 1, then 

(iv) For all {fi,B) E the thermodynamic limit of the susceptibility exists and is given by 


x{fi,B) := lim XN{fi,B) = -^fi{/3,B). 
N^oo OB^ 


(1.13) 


Having investigated the phase diagram of the annealed Ising model on the GRGjv('ir’), we next state 
the SLLN and GLT for the total spin in the following two theorems: 

Theorem 1.2 (Annealed SLLN). Let (Gjv)jv>i ^ sequence 0 /GRGjv('*r’) graphs satisfying Condition 
ll.il then, for all {fi,B) E for any e > 0 there exists a number L = L{e) >0 such that the total spin 
is exponentially concentrated in the form 

for all sufficiently large N, 
where M = M{fi,B) is the annealed magnetization defined in (jl.l2h . 




N 


-M 


> 


Theorem 1.3 (Annealed GLT). Let (Gjv)^>i be a sequence of GRGjv('ur) graphs satisfying Condition 
li.il Then, for all {ffiB) E the total spin satisfies a CLT of the form 


Sn ~ Pn (Sn) V ^ 
y/N 


A^(0,x), 


w.r.t. P, 


iV) 


as N ^ 00 , 


where x = xiP^P) *■5 thermodynamic limit of the annealed susceptibility defined in (11.13^ andj\f{0,a‘^) 
denotes a centered normal random variable with variance cr^. 


The proofs of Theorems 11.11 11.21 and 11.31 all heavily rely on the fact that the annealed GRGjv('ir’) 
gives rise to an inhomogeneous Gurie-Weiss model, which is interesting in its own right. We continue by 
studying the annealed measure on GMjv(2). 
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1.4.2 Results for CMjv(2) 

Our main result for CMjv(2) concerns its thermodynamic limits, a SLLN and a CLT for its total spin, as 
formulated in the following theorems: 

Theorem 1.4 (Thermodynamic limits for the annealed CMjv(2)). Let (G'jv) 7 v>i ® sequence o/CMjv(-2) 
graphs. Then, for all (5 > t), B ^ M, the following hold: 

(i) The annealed pressure exists in the thermodynamic limit N ^ oo and is given by 

if{l3,B) := lim ^^{ld,B) = \og\+{(3,B), 

N^oo 


where 


A+(/3,R) = 


cosh(R) + Y^sinh^(R)”-|-"e^^ 


(a) The magnetization per vertex exists in the limit N ^ oo and is given by 

M{P,B) := lim M^{P,B) = smh(R) ^ 

^ sinh^(R) + 


(1.14) 


Remark 1.1. Since lim 5 ^Q+ M{(d,B) = 0 for all /3 > 0, by definition (jl.Qp we conclude that there is no 
annealed phase transition for CMjv(.2). This is not surprising, since CM^(.2) consists of a collection of 
disjoint cycles, and the Ising model does not have a phase transition in dimension one. 

Next we state the SLLN for the total spin in CMjv(2): 

Theorem 1.5 (Annealed SLLN for CM^(2)). Let (G'jv)Ar>i be a sequence o/CMjv(.2) graphs. Then, for 
all /3 > 0, B G R, for any e > 0 there exists a number L = L(e) > 0 such that the total spin is exponentially 
concentrated in the form 

for all sufficiently large N, 

where M = M(I3,B) is the annealed magnetization defined in (11.141) . 

Finally, we investigate the CLT for CMjv(2): 

Theorem 1.6 (Annealed CLT for CMjv(2)). Let (Gjv) 7 v>i ® sequence o/CMjv(-2) graphs. Then, for 
all P > 0, B ^ M, the total spin satisfies a CLT of the form 






> 


Sn ~ Pn (Pn) v ^ 

y/N ' 


A^(0,x), 


w.r.t. Pff, as N ^ Qo, 


where x = xiP^P) thermodynamic limit of the quenched susceptibility (see Theorem 1.1]) of the 

Ising model on CMjv(-2). Moreover, y(/3,R) is also equal to the susceptibility of the one-dimensional Ising 
model, i.e.. 


xi/3,B)=x'^=\/3,B) 


cosh(R)e 
(sinh(R) + 


Theorems 11.41 [T3] and [L6] are proved in Section |3l Their proofs heavily rely on the fact that CMjv(2) 
consists of a collection of cycles, and the partition function on a cycle can be computed explicitly. 


9 














1.4.3 Results for CMjv(l, 2) 


Our main result for CMjv(l, 2) again concerns its thermodynamic limits, SLLN and CLT for its total spin. 
Some of the quantities involved in the statement of these results are defined in Section [H 


Theorem 1.7 (Thermodynamic limits for the annealed CMjv(l,2)). Let (Gjv) 7 v>i ® sequence of 
CMjv(l, 2) graphs for a given p G (0,1). Then, for all /3 > 0, B £ M, the following hold: 


(i) The annealed pressure exists in the thermodynamic limit N ^ oo and is given by 


:= lim V'iv(/S,-B) 
= logA+(/3,R) + 


1 — p 


logA+{(3,B)+His*,t*), 


(1.15) 


where A^((3,B) is defined in {4-3) below, the function H: [0, x [0,p] —>■ M zs defined in 14-25 ) 

below, and (s*,t*) is the unique maximum point of {s,t) H{s,t) on [0, x [0,p]. 


(ii) The magnetization per vertex exists in the limit N oo, i.e., 


M{(5,B) := hm M^(/?,R) = ^f^if3,B), 

N^oo U-D 


and is given in ^4-46) below. 


(1.16) 


Remark 1.2. (a) Since hmg^Q+ M{/3, B) = 0 for all jd > t) (the explicit expression of the magnetization 
is given in (I4.46P ). we again conclude that there is not an annealed phase transition also for CMjv(l, 2). 
Again this is not surprising, since CMjv(l, 2) consists of a collection of one-dimensional lines and cycles, 
and the one-dimensional Ising model does not have a phase transition. 

(b)W proved a CLT for the number of lines of given lengths in CMjv(l, 2). Leveraging on this result, we 
proved in m thjG Qj'VG'T'CLQG-d (^UGTljChGdj J'OT' ill/G totOiJj S'P't'Tt oj" iJl/G Is'tfhQ TThodGl OTL hOj'VG 

applied the result of to compute also the annealed pressure (Ull), but obtaining a result different form 
(I1.15P . While we are able to see numerically that the two formulas agree, we have no analytic proof that 
they coincide. 


We next state the SLLN for the total spin in CM^(1, 2): 

Theorem 1.8 (Annealed SLLN for CM^(1,2)). Let (G]v)^>i be a sequence o/CMjv(l,-2) graphs. Then, 
for all (3 > 0, R G M, for any e > 0 there exists a number L = L(e) > 0 such that the total spin is 
exponentially concentrated in the form 

e ) < for all sufficiently large N, 




N 


-M 


> 


where M = M(I3,B) is the annealed magnetization defined in (11.161) . 

We finish with the annealed CLT in CMjv(l, 2): 

Theorem 1.9 (Annealed CLT for CMjv(l, 2)). Let (Gjv)jv>i be a sequence o/CMjv(l,^) graphs. Then, 
for all (3 > 0, B £ M, 


Sn ~ Pn ('S'jv) V ^ 
^/N 

where o\ is defined in (f.^. j5[) below. 


Af {0,al) , 


w.r.t. Pjv, as N ^ oo, 


Theorems 11.7111.81 and 11.91 are proved in Section S] 


10 





















1.5 Discussion 

1.5.1 Properties of annealing 

From the results described above, the following general picture emerges on the effect of annealing: 


(i) First of all, in the presence of a ferromagnetic phase transition, annealing can change the critical 

temperature, meaning that We proved this for the rank-1 inhomogeneous graph. For 

the configuration models with vertex degrees at most two that we have analyzed, it holds (3^^^ = 

= oo. We conjecture that in the general case when there is a positive proportion of vertices of 
degree at least 3 and v > 1 (so that there exists a giant component), an annealed positive critical 
temperature exists. We believe that this annealed critical temperature is strictly larger than the 
quenched critical temperature whenever the vertex degrees fluctuate and a positive proportion of 
the vertices have at least degree three. 

(ii) Furthermore, the annealed state satisfies a central limit theorem for the rescaled magnetization, as 
the quenched state does as proved in our previous paper |8]. Unfortunately, we can only prove this 
for certain random graph sequences, but we believe this to be true in general. The variance of the 
annealed CLT and the variance of the quenched CLT are different whenever the degrees are allowed 
to fluctuate. We showed this in the case of the generalized random graph, where they can not be 
ordered because the quenched and annealed critical temperatures are different and the quenched and 
annealed susceptibilities diverge at the critical point. For CMjv(l, 2) having zero critical temperature 
and fluctuating degrees the variances are also different, and we believe the annealed variance to be 
larger than the quenched variance. Unfortunately, we have not been able to prove this. 

(iii) From the analysis of the CMjv(2), we see that both the annealed critical temperature and the 
annealed variance are the same of their quenched counterparts. We conjecture this behavior to 
occur for all random regular graphs. 

(iv) In the GRGjv('if), when the weights (rci)ig[ 7 v] are i.i.d and such such > w) = 
for some r > 1, the annealed partition function satisfies 

(1.17) 


Thus, the effect of annealing of the weights is dramatic, as the pressure becomes infinite for every 
/? > 0. To see (11.171) . we first note that the upper bound is trivial, as H{a) < N(N — l)/2. Thus, it 
suffices to prove a matching lower bound. With Kff the complete graph on N vertices and for a > 0, 


F{GRG^{w) = K^) = 

ij 


>P(u;i G [Af“,2Ar“]Vf G [N]) 

xE[l[pij\wie [iV“,2iV“]Vf G [N] 
ij 


We analyze both terms separately. Firstly, since the weights are i.i.d., 

F{wi G [A^“,2Ar“]Vi G [N]) = P(r(;i G [A^“,2iV“])^ > . 


Secondly, when Wi G 2N°'] for every i, there exists 5 > 0 such that 


Pij > 1 - bN^-\ 
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Therefore 


E 


[llp^j I Wi G [iV“,2iV“]Vi G [AT]] > (l - 57Vi-“) 


N{N-l)/2 


= 1 - 0 ( 1 ), 


when a > 3. Thus, 

Qr,{Z^i(3,B))>Z^^{f3,B)e° 

where Z^^ {(3, B) is the partition function on the complete graph. This proves the claim. 




We will expand on the analysis of the annealed critical behavior of Ising models on generalized ran¬ 
dom graphs in a forthcoming paper [7], where we study critical exponents around the annealed critical 
temperature and we derive non-classical asymptotic laws at criticality. 


1.5.2 CLT proof strategy 

By applying a commonly used strategy |18) . we can prove CLTs for (5'jv)Ar>i by showing that the moment 
generating function of the rescaled total spin Vff = converges in a neighborhood of t = 0 to the 

moment generating function of a centered Gaussian random variable. The convergence can be achieved 
by considering the so-called scaled cumulant generating functions of Sjv, defined as 


Civ(t) = ^ log E [exp (t5jv)], 


(1.18) 


and by proving the convergence of the sequence (c((,(tjv)) 7 V>i for tjv = o(l) to a finite value y, which turns 
out to be the variance of the normal limit. This strategy has been followed in the quenched setting in |8] 
where, specializing E to the relevant measures, the CLT was proved for the Ising model on the whole class 
of locally tree-like random graphs in the random quenched setting, and for the CMjv(2) and CM^(1,2) 
models in the averaged quenched setting. In the former case, the limit c(t) := limAr^oo Cjv(i) can be 
established as a simple consequence of the existence of the random quenched pressure on locally tree-like 
graphs, while the convergence of (c((,(tjv))iV>i follows from the concavity of the first derivatives of the 
cumulant generating functions. In the random quenched setting, this in turn is a consequence of the GHS 
inequality, which holds for the ferromagnetic Boltzmann-Gibbs measure PGm- other hand, under 

the averaged quenched measure this derivative can not be expressed in terms of the averaged quenched 
magnetization to exploit the GHS inequality. Because of that, only the CMjv(2) and CMjv(l,2) setting 
have been treated in j8] explicitly, by exploiting the structure of the graphs and connecting these systems 
to the one-dimensional Ising model. 

A similar scenario is found in this paper, where the approach to the proof of the CLT described above 
is applied to the annealed setting, i.e., with (11.1811 replaced by the annealed cumulant generating function 


Civ(t) 


— logPjv [exp (tSf,)] 


that can be connected to the annealed pressure, since Cjv(t) = '0Ar(/3, B + t) — fjNiP, B), see (II. 8p . 

We will show by an explicit computation that the annealed pressure of GRGjv('IA’) coincides with that of 
an inhomogeneous Curie-Weiss model. From this fact, the thermodynamic limit of the annealed pressure, 
magnetization and susceptibility can be obtained. This again relies on the GHS inequality that is valid 
also for this inhomogeneous ferromagnetic system. Thus, for the generalized random graph, the annealed 
CLT can be proven in a similar way as for the random quenched measure. 


12 






On the other hand, the proofs of the CLT for the configuration models do not follow from the abstract 
argument based on the GHS inequality, since GHS is not available in the general annealed context. Because 
of that, we have to explicitly control the limit (c^(tiv))Ar>i throughout the computation of the annealed 
pressure. It is relatively simple to accomplish this task in the case of the regular CMjv(2) graph consisting 
of cycles only. The fluctuating degree of CMjv(l,2) makes the computation of the pressure and of the 
limit much more involved. CM^(1,2) consists of both lines and cycles. While the cycles give 

a vanishing contribution to the thermodynamic limit, the distribution of the length of the lines has to be 
carefully analyzed and its Gaussian fluctuations appear in the GLT for the total spin. 

1.5.3 Paper organization 

The rest of the paper is organized as follows. In Section [2] we deal with GRGjv('*a’) for which we compute 
the pressure and magnetization in the thermodynamic limit, identify the critical temperature and then 
prove the SLLN and GLT. All of these results rely on the fact that the Ising model on GRGjv('if) in 
the annealed setting turns into an inhomogeneous Gurie-Weiss model. The pressures and GLTs for the 
2-regular configuration model are considered in Section |3] and for the configuration model with vertex 
degrees 1 and 2 in Section 21 In the former case, we show that the variance of the limiting normal variable 
is the susceptibility of the one-dimensional Ising model. In the latter case, which is much more difficult, 
the varying degrees of the vertices affect the pressure and the limiting distribution. In fact, the limiting 
variance is the sum of that of the one-dimensional Ising model and of an extra term emerging from the 
fluctuations of the connected structures of the graph. 

2 Proofs for GRGiv(w) 

GRG In this section, we derive our results for the generalized random graph GRGjv('U^) stated in Theorems 

rmrr^ andrrm 

2.1 Annealed thermodynamic limits: Proof of Theorem 11.11 

The proof is divided into several steps. 


Annealed partition function. We start by analyzing the average of the partition function for GRGjv('a’). 
By remembering that in this random graph the edges are independent and denoting by lij the Bernoulli 
indicator that the edge between vertex i and vertex j is present, we compute 


Qpf (Zjv (/S, B)) = Qjv ( ^ exp [^/3 ^ + B at 

crgQjv i£lN] 






We rewrite 


irSOjv i<j 

Yl + (1 - Pij 

i<j 
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where /3ij and Cij are chosen such that 

e~^Pij + {1 - Pij) = Cije~l^*\ and pij + {1 - pij) = Cije^^K 
Now, by adding and dividing the two equations of the system above, we get 

Cij cosh (Pij) = pij cosh(/3) + (1 - Pij) , ^ log -p\ 

Then, using the symmetry Pij = Pji we arrive at 

Q^{Z^{P,B)) = (n^b) E 

i<j crSOjv 

= G{P)Gi{P) ^ ^ (2.1) 

where G(/3) = Wi^j Cij and Gi(/3) = OieiAr] and we write pa = + wf )- This is the starting 

point of our analysis. We can recognize the r.h.s. as the partition function of an inhomogeneous Ising 
model on the complete graph, where the coupling constant between vertices i and j is equal to Pij. In the 
next step, we analyze this partition function in detail. 


Towards an inhomogeneous Curie-Weiss model. We continue by showing that Pij is close to 
factorizing into a contribution due to i and to j. For this, by a Taylor expansion of x log(l + x), 

Aj = ^log {l+pij{ef^ - 1)) - ^log (l +Pij{e~^ - 1)) 

= - 1 ) - - 1 ) + Oipjj) = s\uh{P)pij + 0{pjj). 

Then, 


Q^{Z^{P,B)) 

= cm Y. 


( 2 . 2 ) 


e[iv] 


(Ti I sinh(/3) Eije[iv] PO<^i<^j+0(Eije[jv] 


(tG^N 


where G 2 {P) = G{P)Gi{P). To control the error in the exponent, we use pij < WiWj/ij.^ and the assump¬ 
tions in Condition II.1[ to obtain 


E p 

i,je[N] 


2 ^ ^ 
ija^aj 


< 


E 

i,j&[N] 


WiWj 


'Eie[7V] 


= o{N). 


Then, 


_ -I . Wj^Wj 

QMiZM{P,B)) = G2iP)e°^^'> Y e-®^“e[iv]<^ie2«“M/3)Ei.j6[A,]— 


iV 


G2(/3)e°(^) Y e^^*e[^T'e2^ 


1 sinh(,S) /v-' 


When Wi = w for all i, so that GRGjv('u;) is the Erdos-Renyi random graph, we retrieve the Curie-Weiss 
model at inverse temperature P' = sinh(/3)u;. In our inhomogeneous setting, we obtain an inhomogeneous 
Curie-Weiss model that we will analyze next. 
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Analysis of the inhomogeneous Curie-Weiss model. We use the Hubbard-Stratonovich identity, 
i.e., we write e* = E[e*^], with Z standard Gaussian. Then, we find 


f7Gr2jv 

N 

= G2(/3)e°(^)2^EjjJcosh 




sinh (/3) 


2=1 


= G2(/3)e'’(^)2^E 


(-N 

exp I log cosh 


N 


2=1 


WiZ + B 


I sinh(/3) 


WiZ + B 


)} 


We rewrite the sum in the exponential, using the fact that = wij^, where we recall that Ijv is a uniform 
vertex in [N], to obtain 


Q^{Z^{/3,B))=G2{/3)e°^^h^E 
= G2{/3)e°^^h^E 


|aE log cosh 


exp 




sinh(/?) 

AE[Wjv] 


Wj,Z + B 


M. 


where 


Fff{z) = E 


log cosh 


I sinh (/?) 


Wt,z + B 


, E[W^] 

Here we emphasize the fact that in (|2.3I) . the expectation is w.r.t. Wj^ only. 

We continue by analyzing Fr^{z). We claim that, uniformly for \z\ < a and any a < oo. 


sup \F^{z) - F{z)\ = o(l), 

|2:|<a 


(2.3) 


(2.4) 


where 


F{z) = E log cosh 


I sinh (/3) 
E[W] 


Wz + B 


To see (12.41) . we note that Fj^{z) —>■ F{z) for every z fixed by Condition ll.ll ai-fbb and the fact that 
logcosh(x) < |x|. Further, 


sinh(/3) 

\Fj^[z)\ < . . E 


E[W^ 


tanh 


I sinh (/?) 
E[1T^] 


W^z + b]w^ 


< sinh(/3). 


since tanh(x) < 1 for all x, so that \F'^(z)\ is uniformly bounded in N and z. Therefore, (T)v)Af>i 
forms a uniformly equicontinuous family of functions, so that (12.4p follows from Arzela-Ascoli. Since 
Fn{z) < sinh(/3)|z|, it further follows that, for a > 4sinh(/3), 


E 






{|Z|>aViV} 


< E 

= 2E 
2 


jvWsinh(^)|Z| j 
^\/N sinh(/3)Z j 


{\Z\>aVN} 


{Z>aVN} 


sinh(l3)z^-z'^/2^^ 


2 roo 

JaVN 

poo 

^ ^asinh{/3)N-a‘^N/2 / ^VN{sinh{l3)-a)x ^ ^-a'^N/A 

Jo 
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which, for a sufficiently large, is negligible compared to E 



^{\Z\<aVN} 


We conclude that 




G2(/?)e°(^)2^E 


■ NF 
e 



( 1 + 0 ( 1 )). 


(2.5) 


A large deviation analysis. The expectation in (j2.5p is an expectation of an exponential functional, 
to which we apply large deviation machinery. The Gaussian variable Z/y/N satisfies a large deviation 

principle with rate function I{z) = 2^/2 and speed iV, because Z/y/N = ^ {Zi + ... + Z^^), where (^i)ig[Ar] 
are i.i.d. standard Gaussian variables. Using Varadhan’s Lemma and the fact that z ^ F{z) is continuous, 
we calculate the thermodynamic limit of the pressure as 


lim ^ log Qjv {Zn (/?, B)) = log 2 + lim ^ log G 2 (^) + sup [F{z) 

N—^oo N^oo 2 


- m 


= log 2 + a (/3) 


+ sup 

Z 


E 


log cosh 


I sinh(/3) 
E[1U] 


Wz + B 


-1 
2 . 


( 2 . 6 ) 


where a (/3) = limv^oo ^ logG 2 (/3). The equation that defines the supremum is 


z* = z*{j5,B) = E tanh 


/ smh(/3) 

EIW] 


+ B 


sinh (p) 
E[W] 


W 


(2.7) 


and the annealed pressure is obtained by substituting the supremum point z* in the right hand side of 
(|2.6p as 


tp(/3, B) = log 2 + a(/3) + E log cosh 
This completes the proof of Theorem 0(1). 


/sinh (/3), 
E[1U] 


Wz*{P,B) + b] -z*{p,Bfl2. 


( 2 . 8 ) 


The critical inverse temperature. To identify /3™ as stated in Theorem ll.ll iiL we evaluate (12.7p 
when B '\0 to obtain 


z* = H{z*) where 


H{z) = E 


tanh 


/ / sinh (/3) 
E[W] 



/sinh (/3) 
E[W] 


W 


(2.9) 


We investigate the solutions of z* = H{z*) in (j2.9p . We note that z i—)• H(z) is an increasing and 
concave function in [0,oo). When H'{0) > 1, we have three solutions of (12.91) . i.e., Fz* and 0, where 
2 * = z*(/3,0"^) > 0. When H'{0) < 1, instead, z* = 0 is the only solution. This leads us to compute that 

H (0) = sinh {(5) = sinh (/3) v. 

Thus, the annealed critical temperature /?“ satisfies sinh(/?“) = X/v. Since tanh (/3c'^) = l/i/, and 
tanh(x) < sinh(x) Vx > 0, we obtain /3c'^ > (3^^, unless when u = oo, in which case /3“ = /3c'^ = 0. 
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Thermodynamic limit of the magnetization. To prove the existence of the magnetization in the 
thermodynamic limit stated in Theorem ll.ll iih we follow the strategy used in |14j . We use the following 
lemma: 


Lemma 2.1. Let {fn)n>i be a sequence of functions that are twice differentiable in x. Assume that 

(a) lim^^oo fn{x) = f{x) for some function y i—>■ f(y) that is differentiable in x; 

(b) -§^fn{x) is monotone in [x — h, x + h] for all n > 1 and some h > 0. 

Then, 

lim ^fnix) = ^fix)- 
n^oo dX dX 

We apply Lemma [2.II with n = N and equal to B V’iv(/3) B). We verify the conditions in Lemma [2.11 
and start by noting that 

Mj^iffB) = Pj,{Sj,/N) = ^f^{/3,B), 

and limjv^oo V'iv(/3,1?) = ipiffB) by Theorem II.If ij with B i—>■ Mff{l3,B) non-decreasing: 


d 


dB 

Thus, we can indeed conclude that 


— B) = - [P^ (52 ) - P^ {S^f 


> 0 . 


MiffB) = lim Mr^iffB) = lim —f^iffB) = —f{ffB). 

N—foo N^oo uJd uJd 

The limit magnetization M(l3,B) can be explicitly computed by taking the derivative of if{ffB), 
and using the hxed point equation (|2.7p . to obtain 


M{ffB) = E 


tanh 


I sinh (/?) 
E[W] 


Wz* + B 


Thermodynamic limit of the susceptibility. Finally, the thermodynamic limit of the susceptibility 
in Theorem ll.ll iv) is proved using Lemma 12.11 by combining Theorem ll.ll iii and the fact that B 1 —)• 
B) is non-increasing by the GHS inequality. Indeed, by the explicit computation in (j2.ip . we 
see that the annealed partition function can be viewed as the partition function of an inhomogeneous 
Curie-Weiss model, where the held is homogeneous and the coupling constants depend on the edges. Since 
such an inhomogeneous Ising model also satishes the GHS inequality, the same follows for the annealed 
partition function for GRG^(rc). Therefore, 

Q^AlMiffB) = 


□ 
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2.2 Annealed SLLN and CLT: Proofs of Theorems [TT^ and [T73] 

With Theorem 11.11 in hand, we now have all the hypotheses to prove Theorems 11.21 and 11.31 following the 
strategy used for the random quenched setting in Sections 2.2 and 2.3 of |8] verbatim. Indeed, for the 
proof of the annealed SLLN, referring to [H Section 2.2], we obtain the existence of the thermodynamic 
limit of the annealed cumulant generating function 

Civ(t) = ^logPv [exp(tS'iv)] = ipN{P,B + t) - 

by Theorem II.H i). Then, from |18l Theorem 11.6.3] and Theorem I l.ll iii we conclude the proof. 

To prove the annealed CLT (see [U Section 2.3] for the proof of the random quenched CLT) we need 
the existence in the thermodynamic limit of pressure, magnetization and susceptibility given by Theorem 
11. II together with the CHS inequality that is still true in the annealed setting thanks to the mapping to 
the inhomogeneous Curie-Weiss model. 


3 Proofs for CM^(2) 

In this section we prove the CLT with respect to the annealed measure for the 2-regular random graph. 
We start by computing the annealed pressure using the partition functions for the one-dimensional Ising 
model with periodic boundary conditions. 


3.1 Annealed thermodynamic limits and SLLN: Proofs of Theorems fTT^ and fT75] 

From our previous paper [8], we remember that any 2-regular random graph is formed by cycles only. 
Thus, as in |8], denoting the random number of cycles in the graph by we can enumerate them in an 
arbitrary order from 1 to and call the length (i.e., the number of vertices) of the ith cycle. The 

random variable has distribution given by 

N 

Ki = Y,Ij, (3.1) 

i=i 

where Ij are independent Bernoulli variables given by 

Ij = Bern (- | . (3.2) 

^ \2N-2j + lJ ^ ’ 

See [8] for a proof of this fact. Since the random graph splits into (disjoint) cycles, its quenched partition 
function factorizes into the product of the partition functions of each cycle. Therefore, 

Z,(/3,B) = n^Li(q(/3,^). (3.3) 

i=\ 

By [g Section 3.1], we have that the partition function of the one-dimensional Ising model with periodic 
boundary conditions is given by 

Z^‘H/3,11) = A(^'(/3,A) + A(!(/3,A), (3.4) 


where 


A±(/3,S) = e^ 


cosh(A) ± y^sinh^^Ay^t-"e^^ 


(3.5) 
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so we can write 


Because /3 > 0, we have 0 < A_(/3, B) < A+(/3, B), so that, for every i, 

As a result, we can bound the the pressure as follows: 

T^t T^t T^t 

J^N 

i=l i=l i=l 

and, since XliJi ^^{i) = N, we finally obtain 

A^(/3,i?) < Z^{P,B) < 2^nxI{P,B). (3.6) 

The thermodynamic limit of the annealed pressure 'ipniP, B), defined in (|1.8p . can be computed along the 
same lines of the averaged quenched one in [8]. Indeed, by applying the monotone operator N~^ log(Q]v(‘)) 
to (j3.6l) and using the fact that A+(/3, i?) is non random, we obtain 

logA+(^,5) < ^n{P,B) < ^log (^Qjv +logA+(^,5). 


Now using the fact 


1 

iv 


log (Qn ( 2 ^^)) 


1 ^ 

^logllQ- (2^0 


N 

1 

iV 


2 = 1 
N 

log]J 


2=1 


N 




2N -2i +I 
1 


+ 1 - 


1 


2N -2i +I 


2 = 1 


2A^ - 2i + 1 


N^oo 


0, 


(3.7) 


we conclude that the annealed pressure of CMjv(2) coincides with the pressure of the one-dimensional 
Ising model {P, B), i.e., 

P{P,B) = P^^=\P,B) = \ogX+{P,B). (3.8) 

Moreover, it also agrees with the averaged and random quenched pressures [8], i.e., 


ip{P, B) = ip{P, B) = pjiP, B), 


where 


'p{P,B) := lim ^Q^{ln Zn{P, B)) 

N^oo 


and 


P’W,B) := lim ^lnZAr(^,B) . 
N-^oo 


It also straightforwardly follows that the annealed cumulant generating function of CM^(2) coincides with 
the random and averaged quenched ones |8] i.e.. 


c{t) = c{t) = c{t) = logX+{P,B + t) -logX+{P,B). 


(3.9) 
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The existence of the magnetization in the thermodynamic limit lTheorem ll.4l fiili can be proved, as in the 
previous section, using Lemma l2.1l and the existence of the thermodynamic limit of the pressure ()3.8p . so 
we obtain 




A 




sinh(i?) 

yAih^(fi)”T"e^^ 


as required. 


□ 


Proof of Theorem L.5i ' The proof, as for Theorem 11.21 follows immediately from the existence of the 
annealed pressure in the thermodynamic limit and its differentiability with respect to B. See also |8l 
Section 2.2]. □ 


3.2 Annealed CLT: Proof of Theorem 11.61 

To prove the CLT in the annealed setting, we follow the strategy used in jS] for the averaged quenched CLT. 


Rewrite in terms of cumulant generating functions. Using the annealed cumulant generating 
function and using a Taylor expansion, we write 


logPjv 


exp 



■ 7 ^ ) 



(3.10) 


where G [0,t/y/N]. Then the aim is to prove that limAr^oo ^(^iv) exists as a finite limit. 

By expressing Cjv(t) in terms of and using (13.91) . we can compute the difference as 


c^(t) - c(t) = — log 


Qi 


I (A+(/3,B+t))" 


Qi 


I (A+(/3,B))" 


1 




" <3»(n3(i + (’B)‘»<-')) _ 


where, as in [8], we have defined 


rg = r{l3,B) = 


A-(/3,R) 

A+(/3,R) • 


Then 


Cjv(t) = logA+(/?,R + t) -logA+(/?,R) + 


Qn ^ 

nS(i + (^B+t)"-«) 

)i 

Qn 

(n*7L (1 + 



(3.11) 


(3.12) 


Our aim is to show that the double derivative arises from the first term only, the second derivative of the 
last term vanishes. 


Computation of the second derivative of the cumulant generating function. The second deriva¬ 
tive of (|3.12l) is 


= ^logA+(/3,R+ t) + 


1 


ND^it) 


lN(t) + + 


Dpf(f) 


(3.13) 
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where 






2=1 


J^N 




i=i 

j¥=i 


^N ^N 


^N 


(l + (?’s+t)'^^^*^) 


i=l j=l 

j¥^i 

ffi 


1=1 

i¥=hj 




Ill^it) = J] (1 + )]' 


2=1 


j = l 
j¥=i 


J^N 




2=1 


Uniform bound of the averaged normalized partition function. To analyze the contributions 
above we show that the averaged normalized partition function of CMAr( 2 ) is uniformly bounded: 

Lemma 3.1 (The partition function on tori). For every 7 < 1 and a € (0,oo), there exists a constant 
A = ^ 4 ( 0 , 7 ) such that, uniformly in N, 


K 


t 

N 




2=1 


Proof. Denote = QN\\\f=i + 07 ^^^*^)]- For the proof we use induction in N. The induction 
hypothesis is that there exists an j4 > 1 such that 

Fix M >1 large. We note that we can hx A so large that the inequality is trivially satisfied for N < M. 
To advance the induction hypothesis we first derive a recursion relation for Zi\j. We have 


Zn 


N Kf, 

^ g^(Ljv(i) = 0 Qv(n (1 + 


1=1 


2 = 1 


Ln{ 1) = l) 


N 

J2Qn{Ln{1)=1) (l + aV)^iv-, 

1=1 


(3.15) 


Indeed, the average of HiJi conditioned on Tjv(l) = I, reduces to the average on a CMjv(2) 

graph with N — l vertices of a similar product. This average gives rise to the factor 2^jv-i in (|3.15p . while the 
term corresponding to the hrst cycle is factorized, being (l + 07 *). Substituting the induction hypothesis 
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into (13.151) leads to 


< 


It is not hard to see that 


Z^<A^Q„{L^{l)=l){l+ aj‘) (1 - 

N N 

Aj2QNiLN{l) = 0(1 + «V) -AJ2QALn{ 1) = l)-j^== 

1=1 1=1 ^ 1 


N 


Y,Q^{Ln{ 1) = IW <c/{N+ 1), 


1=1 

while there exists a constant 9 > 1 such that 

N 


Y,Qn{Ln{1) = 1) 


> 


9 


1=1 


-I+ i ~ ^fWTi 


(3.16) 


Indeed, by m Exercise 4.1], or an explicit computation, LAr(l)/iV —^ T, where T has density frix) 
given by 


frix) = 


1 


2y/T^' 


Therefore, rewriting the sum in (I3.16P we have: 
N 


(3iv(.b7v(i) = O' 


1 


1 


-Qi' 


1 


W^TTT WTT L^i-LAr(i)/(iV + i)J’ 

and by Fatou’s Lemma and weak convergence, we obtain 


lim inf Qpf 

N^OO 


r 1 

> E 

r 1 1 

^^l-LAr(l)/(iV + l)J 




> 1 . 


Since we can assume that N > M, which is sufficiently large, we thus obtain (I3.16p . Thus, 

1 


^ ^1 1 + 


A^ + 1 2^iV + 1 


< ^ 1 - 


2^N + 1 


when N is sufficiently large. This advances the induction hypothesis and completes the proof of the 
lemma. □ 


Analysis of the second derivative of the cumulant generating function. Armed with Lemma 
13.1[ it is now easy to show that all the contributions in the second term of the r.h.s. of (13.131) indeed 
vanish on a sequence = o(l). To see this, let t > 0 and (tjv)Af>i a sequence of real numbers such that 
tjv £ [0,t/vOV]. We consider first the term /jv(tjv). As in (8] Lemma 3.1], there exists a constant C > 0 
such that 


J^N 


\LAi){LAi) - < C • 


2=1 
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since < 1. Then, using the Cauchy-Schwarz inequality and recalling that 0 < r < 1, we obtain 

i=l 

< c-Qn(^{kIj) ^ • Qiv(^ JJ (i + (?'s+tjv)^^^“') ) 

i=l 

< CQ^ ( ( n (1 + 

i=l 

Using Lemma l3.II with a = 3 and 7 = we conclude that 

i=l 

Finally, since > 1, 

|^iv(tjv)| ^ C • • log Nj^ Q 

ND^itr,) ~ N ^ • 

Similar computations allow us to estimate /Ijv(iiv) and ///jv(tjv) to obtain 


lim 

TV—>-00 


//jv(tAf) 

N Z)jv(tjv) 


= 0 , 


lim 

TV—>-00 


IIlN{tN) 


= 0 . 


Completion of the proof of Theorem 11.61 Having proved that 


1 


lim 

N^oo ND!^{tf 


-^iv(tjv) + -^-^iv(tAr) + 


m n(J^n) 


= 0 


A*iv(tjv) 

the combination of (|3.10p and (13.131) yields the proof of the annealed CLT, i.e., 


lim logPjv 

N^OO 


exp t 


<Sjv — Pn{Sn) 

^/N 




■log \+{/3,B + t) 


2 dP 

cosh(H)e“^^ 
2 (sinh(i?) + 


t =0 


Therefore, we conclude that the annealed CLT has the same variance as in averaged quenched case [ 8 ], 
i.e., the variance in both cases is the susceptibility of the one-dimensional Ising model. □ 


4 Proofs for CMjv(l, 2) 

In this section, we consider the Configuration Model CMjv(l, 2), introduced in Section 1.2. In this graph, 
the connected components are either cycles or tori (which we indicate by a superscript (f)) connecting 
vertices of degree 2 , or lines (indicated by a superscript {1)) having vertices of degree 2 between two vertices 
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of degree 1. In order to state some properties of the number of lines and tori, we need to introduce some 
notation. By taking p € (0,1), let us define the number of vertex of degree 1 and 2 by 

ni := #{i E [N]:di = 1} = N - [pN\, := # {z E [N]: di = 2} = [pN\, 

and the total degree of the graph by 

= 2^2+ ni = N+[pN\. (4.1) 

ielAf] 

Then, the number of edges is given by /2. Let us also denote by Kff the number of connected components 
in the graph and by and the number lines and tori. Obviously, 

I ■ 

Because every line uses up two vertices of degree 1, the number of lines is given by ni/2, i.e., K^'' = 
{N — [pA^J)/2 a.s.. Regarding the number of cycles, we have that has the same distribution of 
where N is the (random) number of vertices with degree 2 that do not belong to any line and is the 
number of tori on this set of vertices. Then, since this subset forms a CMjv(2) graph, we can apply [HI 
(3.16) in Section 3.2], obtaining that /N —0, so that also 

K^/N -^{l-p)/2. 


Denoting the length (i.e. the number of vertices) in the zth line and jth torus (for an arbitrary labeling) 
by and the partition function can be computed as 


K 


(0 


K 


(t) 


ZAfl.B) = n ■ n 


2 = 1 


2=1 




(4.2) 


where, by (13.4h . 


Z 


(t) (a 


while the partition function on each line is obtained using the partition function on one-dimensional Ising 
model with free boundary condition [HI Section 3.1] as 


where 


A± = A±i/3,B) = 


=A+\^ +A. 

e-2/3g±s ^ ± 2e-^{\+ - e^+^) 

[e-2/3 + (A+ - e/3+^)2]A± ' 

This is the starting point of our analysis of the annealed Ising measure on CMjv(l, 2). 

4.1 Annealed CLT: proof of Theorem 11.91 

In order to prove the CLT in the annealed setting, we will show that 


(4.3) 


lim Pjv 

N^OO 


exp - Pn{S^ 


= exp(cr2t^/2), t E M. 


(4.4) 
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From now on, to alleviate notation we will omit the dependence on (3 and abbreviate = B + We 
start by writing 


Pn 


exp 



QAZ^{B)] Q^[^NFs{p(^y)+NE^{B^ 


where, by |8], 


Fb{p^^^) = log A+(5) + Y.Pr' log (^+(^) + A-{B) {r{B)y) , 

l>2 


Er,{B) 


1 

iv 


^log (l + r{B)^N^^py 


i=l 


(4.5) 


(4.6) 


with r{B) 
by 


rs defined in (|3.1ip and = 



the empirical distribution of the lines lengths given 



1 

N 


K 


(0 

N 


2=1 


(4.7) 


Analysis of the annealed partition function. We have 

CX) ^ 

eiVF,(p(-)) ^ (A+(S))^ J] (yl+(i?) + yi_(s) (riB)y) ‘ 

1=2 

OO OO 

={\+{B)f n n (i+«(^) ' • (4.8) 

1=2 1=2 


where 


and Ni 


a{B) 


A-{B) 

A+{B) 


Np^j^'^ is the number of lines of length 1. We rewrite the second factor in (14.81) as 


ll{A+{B)f^ = {A+{B)rB2 ^ 
1=2 


since X )/>2 ~ ni/2. Therefore, we arrive at 


OO OO 

gAfFB(p(^)) ^ A(^'(S)2l+^/^(5) J| (^l + a(B)r'(S)) ' = A(^^(5)^+^/^(5) JJq(5)^', 

1=2 1=2 

where we define 

ci{B) ■.= l + a{B)r\B) . 

Next, define 

M^ = N -^INi 

l>2 


(4.9) 
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for the number of vertices that are not part of a line. Then, denoting by (B) the partition function of 

CM^(2), 

Q^iZ^m Q,[e^Fsip(-))z^^liB)] ^ ^ 

X^{B^)Al^/\B^)Q4Z^^'>jB^)YlZ2Ci{B^)'^q 
X^{B)Al^^\B)Q4Z<;^liB)UZ2Ci{Br^] 

where we write 

Z^^\B) = X^^Z^^\B). 

Asymptotic behavior of the annealed partition function. The key result for the proof of Theorem 
[TQ] is the following proposition that establishes the exponential growth of the annealed partition function 
with polynomial corrections: 

Proposition 4.1. The following holds true: 

(a) For B ^ 0, there exist I = I{B) and J = J{B) such that, as N ^ oo, 

CO 

Q^[Z^^l{B)\{ci{Bf^]=J{B)eBB'^^{l + o{l)). (4.11) 

1=2 

The function B J{B) is continuous, while B I(B) is infinitely differentiable. 

(b) Given t G M there exist I = I(t) and J such that, as N ^ oo, 

n <=‘ {■^) 1 = + »(1)). (4.12) 

The function 1 1 —)• 7(t) is infinitely differentiable. 


Proof of Theorem \l.y\ subject to Proposition ^. 1\ We start proving the theorem for B 0. We substitute 
()4.1ip into (I4.10p to arrive at 


Qn[Zn{Bn)] 

QAZ^B)] 


( 1 + 0 ( 1 )) 
(1 + 0 ( 1 )) 




X^{B)Al^^^{B)J{B)eBB)N 

ni/2 


/ A+(lljv)' 


V K{B) . 

/ V A+{B) ) 


NiI{BN)-l{B)) 


(4.13) 


where we use the fact that B i—>• J{B) is continuous to obtain that J{Bpf) = (1 + o(l))J(i?). We can next 
use the differentiability of B I(B) and the fact that Hjv = B F t/y/N to expand out 


QN[ZN{Br,)\ . 4VAr[|.logA+(i?+t)|t^o+|7f logA+(B+t)|t=o+^/(S+dlt=o] 

Qr,[Zr,m + 

X (4.14) 
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where 


^2 — log A+(i? + t)|t=o +-^—+ t)\t=o +^^I{B+ t)\t=Q. (4-15) 

Since 

h{S^) =N[^^\og\+{B + t)\t=o + ^^^\og A+{B + t)\t=o + + t)|t=o] 

+ o(Viv), 

then (I4.14jl implies (14.4^ . thus proving the theorem in the case B ^ 0. 

For i? = 0, in a similar way now using (|4.12p . we get 

QN[ZN{t/VN)] _ , iy7v[^logA+(t)|t=o+|i^logA+p)|t=o+^/p)|t=o] 

where 

^2 = ^logA+(t)|t=o+ ^ ^^ ^log^+(t)|f=o + ^/(t)|f=o- 

□ 

Strategy to prove asymptotic behavior. The remainder of this section is devoted to the proof of 
Proposition 14.11 We use the law of total probability to write 

CO 

Q4Z^^liB)llci{B)^^] 

1=2 

n 2 oo 

= ^E^[Z^\B)]Q4llci{B)^^ I M^ = m]Q^{M^ = m), (4.16) 

m =0 1=2 

where we denote by the symbol the expectation with respect to an independent CMm(2). 

Our aim is to prove that the asymptotic behavior of (14.161) is essentially dominated by the term with 
m = 0, which gives the exponential growth stated in Proposition 14. II To achieve a full control 

we analyze in the following the three contributions whose product gives rise to the summand of (I4.16p : 

i) Em[Zm\B)]: this is subdominant in the limit N ^ oo since, by Lemma l3.ll sup^Em[Zm (B)] is 
bounded. Therefore it will appear only in the prefactor J{B). 

ii) = m): we study the distribution of the number of vertices in tori Mjv in Lemma l4.lt in 
particular we prove the existence of a limiting distribution function in the limit N ^ oo. 

iii) Qv [ 01^2 I = ITT-]- this is rewritten explicitly in Lemma 14.21 and its asymptotics is 
computed in Lemmata 14.31 and 14.41 


The number of vertices in tori. We start by analyzing the random variable Mjv representing the 
number of vertices belonging to tori. 
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Lemma 4.1 (The number of vertices in tori). When N — oo, there exists a random variable M such that 




■D 


M. 


Further, 


Qn{M!„ = m) 

1 


(ni + 2n2 - 1)!! 


("') 

\m J 


2>^2 ™’(7),2 — m)!(ni — l)!!(2m — 1)!! 

2m 


ni/2 + n 2 — m — 1 
n 2 — m 


^ 2n2 (ni - l)!!n2! f /ni/2 + n2 - m - 1 

(ni + 2n2 — 1 )!! \mj\ n2 — m 

Proof. Number the vertices of degree 2 in an arbitrary way. We write 




n2 


where 




/=i 


2=1 


(4.17) 


and Ji{l) is the indicator that vertex z is in a cycle of length I of which vertex i has the smallest label. We 
compute that 

QnWi^^] = ^(3]v(vertex 1 is in cycle of length 1) —^ ^(2p/(l +p))^ = A;. 

It is not hard to see, along the lines of m Proposition 7.12], that )/>i converges in distribution 
to a collection of independent Poisson random variables {Pi)i>i with parameters (A;)/>i. Further, since 
< ^(n 2 /^]v)^ which decays exponentially, the contribution from large I equals zero whp, i.e., 
Qn(^1 > T: > 0) is small uniformly in N for T large. This shows that 


Mjv -^'^lPi = M. 

l>i 


Note that 

OO OD 

Q^[b^l] =Y[Q^[b^Pt] = (4.18) 

1=1 l=i 

which is finite only when 6 < (1 + p)/{2p). 

To prove (|4.17p . we note that 

Qn{Mm =m)= 2 n 2 - 

where N{ni,n 2 ,m) is the number of ways in which the half-edges can be paired such that there are 
precisely m degree 2 vertices in cycles. We claim that 

N{ni,n 2 , m) = 2 ^^-^{n 2 - m)!(ni - l)!!(2m - 1)!! + ^^2 - m - 1A 

\m J \ n 2 — m J 

For this, note that 

(1) there are (()^) ways to choose the m vertices of degree 2 that are in cycles; 

(2) there are (2m — 1)!! ways to pair the half-edges that are incident to vertices in cycles; 
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(3) there are (ni — 1)!! ways to pair the vertices of degree 1 (and this corresponds to the pairing of 
degree 1 vertices in lines); 

(4) there are (n 2 — m)\ ways to order the vertices that are in lines; 

(5) there are 2 ways to attach the half-edges of a degree 2 vertex inside a line, and there are in total 

n 2 — m degree 2 vertices in lines, giving ways to attach their half-edges; and 

(6) finally, there are ways to create ni/2 lines with n 2 — m vertices of degree 2. 

Multiplying these numbers out gives (I4.20p . This completes the proof of Lemma l4.ll □ 


Combinatorial expression of the partition function. To perform the asymptotic analysis of the 
partition function I = m], we rewrite it as double sum in Lemma 14.21 and then we 

investigate the asymptotics of the summand in Lemma 14.31 by Stirling’s formula. Finally, in Lemma |4.4[ 
we use the Laplace method to estimate the asymptotics of the double sum. 


Lemma 4.2 (Generating function of number of lines in CMjv(l,2)). For every a,r, for ci 
every I >2, 

oo »il/2n2—m 

Qn[Wci{B)^'^ \ = m] = ^ X] - ■m-), 

1=2 £=0 k=0 


where 


B 


il\n2-m) 



/ni/2—i-\-n2—m—k—l\ 
k \ k ) \ n2-m-k ) 

^ ^ /ni/2-\-n2—'in—l\ 


1 -|- ar^ for 


(4.21) 


(4.22) 


Proof. When = m, we have that n 2 — m vertices of degree 2 have to be divided over ni/2 lines. 
Number the lines as 1,..., ni/2 in an arbitrary way. Denote the number of degree 2 vertices in line j by 
Yj and rewrite 


Qn[]Jci{B)^‘ I M. 


1=2 


ni/2 

N =m] = ^ Qjv (Tl = ii,..., = ini/ 2 ) JJ 

(^Iv* 72) j ^ 


where (ii,... ^ini/ 2 ) is such that n + • • • + ini /2 = ^2 — m. Let [^i/2] ,ni/2}, and expand out 

to obtain 


giv(u = *i,...,h;,/2 = v/ 2 ) E 

(*i,...,*nj/2) rc[ni/2] j'er 


(4.23) 


where the sum over T is over all subsets of [ni/2]. We denote 


ni /2 


N, 


ni,n2—m 


= #{(ii, ... ,ini/2): ij > 0 Vj and ^ ij = n 2 - m} 

i=i 

/ni/2 -|- 772 — m — 

\ n 2 — m 
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so that (14.231) is equal to 


E 


n2—m 


(Ov4ni/2) '' 


ni/2+n2—m— 


E 


1 


rC[ni/2] k=0 

nil2 


fnil2+n2—m— 


(n,."4ni/2) '' n2-m 

1 

fni/2+n2—m—l\ 

\ 712—ra / 

ni /2 
X 


lyE 


ni/2 

£ 


712—771 

{ar^Y E ^’"hiii+-+ie=k)} 

k=0 




f=0 

ni/2 n2—m 


k=0 

2\f.k 


— m — k — 1 
— m — k 


—/ /n\ 

E E 


^=0 /c=0 


/£+A:—1\ /7ii/2—i-\-n2—7n—k—l\ 

\ k ) \ 772 — 771—k / 

/ni/2+n2—m—1\ 

\ 712-777 ) 


□ 


Asymptotics by Stirling’s formula. We continue the analysis by investigating the asymptotics of 
BYY{n 2 ) in (14.221) when £, k and n 2 are of the same asymptotic order. To alleviate the notation we write 
^[a\ [b] when a, b are not necessarily integers. 

Lemma 4.3 (Asymptotics of i?£^\n 2 )). Let Dp = [0, (1 — p)/2] x [0,p]. For external fields B 0, there 
exists a function H{s,t) continuous in Dp and smooth in Dp (the interior of Dp) and a function C{s,t) 
smooth in Dp, such that 

^^v{NH{s,t)} {1 +o{l)), as N ^ oo, (4.24) 

Moreover, H{s,t) is strictly concave on its domain Dp and its (unique) maximum point {s* ,t*) lies in the 
interior Dp. In Dp, the functions are defined as follows: 


H{s,t) 


= {l-p) log “ 2s log(s) - 2 log 

+ s log(ar^) + t log(r) + (s + t) log(s + t) - t log(t) 

- {p-t)log{p-t) - log +plog(p), 


and 

C(s,t) 

Finally, uniformly in {s,t) G Dp, 


1 ^Y-^-s-t){s + t~p 


< CN^/^e^v{NH{s,t)} . 


(4.25) 


(4.26) 


(4.27) 
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Proof. Using Stirling’s approximation in the form n! = e ”re"'\/27rn (1 + o(l)) for n large, taking a, 6 E N, 
we can rewrite the binomial coefficients as 



^[6 log h—a log a—{h—a) log (b—a)]n 


Vb{l + o(l)) 

y/ay/b — ay/2Tm 


Plugging the previous formula into (14.221) . then (I4.24p follows. By inspection, H{s,t) and C{s,t) are 
smooth functions in D° for B ^ 0. The function H{s,t) can be further extended by continuity to the 
boundary dDp of Dp, while C(s,t) cannot be defined in Dp \ D° since it is unbounded there. In order to 
prove concavity of H(s,t), we check that its Hessian matrix Q(s,t) is negative definite on each point of 
Dp. For this, we compute the Hessian Q(s,t) as 

/ 1 I 1 _?_2 1 , 1 \ 

^-S+p-t s + t s 1^-s+p-t S + t 

1 1 1 1 1 1 • 

\ ^-s+p-t~^ s + t - s + p- S+ t p-t 1 j 

The eigenvalues pL~^ and of Q{s,t) are 




+ 


k^-s-t s + t ^-s s 


p — t t 


± 




2n 


p — t t 


1+P 


-s-t s+t 


We can easily see that fi~ < 0, and in order to show that also is negative, we observe that the 
determinant of the Hessian matrix is positive. Therefore, H{s,t) is strictly concave in D° and, by conti¬ 
nuity, concave in Dp. Concavity implies that {s,t) H{s,t) has a unique global maximum in Dp, the 
uniqueness follows by strict concavity and the fact that the maximizer is not on the boundary. In order 
to hnd (s*,t*) := aigmaxH(s,t), and to prove that it lies in Dp, we calculate 


dH{s, t) 
ds 


dH{s, t) 

di 


-I- log(ar^), 
log(s + t)- log(f) 


log -s-t^ + log(s + t) - 2 log(s) 

/1 -|- p \ 

log ( — -S - t j + log(p -t)+ log(r), 


so that {s*,t*) is a solution of the system 




(^7^) (^+b 

{p-t) (s+t) _ 1 

* (^-"-7 


(4.28) 

(4.29) 


(4.30) 


Since B ^ 0, and then both ar^ and r are finite and larger than zero, it easy to see from (|4.30p that the 
maximum point cannot be attained on the boundary. 

The proof of (I4.27P follows similarly, now using that e“”'n”\/27rn < n! < 
e“"'n”\/2'7rn(l + j^^) for every n > 1. The power of N is needed to make the estimate uniform, e.g., by 
bounding s((l — p)/2 — s) > c/N uniformly for s > \/N. □ 
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Asymptotics by Laplace method. In the next Lemma we compnte the asymptotic behavior of (I4.2ip 
by using the discrete analogue of Laplace method. 

Lemma 4.4 (Asymptotics of (5jv[nz^2 \ Mj^ = m]). For every m > 0 fixed, 

CO 

\M^ = m] (4.31) 

1=2 

= 27r (det Q(s*, ))"5 (ft*)’" e^p{NH{s* + o(l)), 

where {s,t) e-)• H(s,t) and {s,t) i—>■ C{s,t) are defined in Lemma \4-d\ is the maximum point of 

H(s,t), Q(s,t) is the Hessian matrix of H and 



Proof. We start by proving (|4.3ip for m = 0. Due to Lemma [4.31 we may estimate the asymptotic behavior 
of the double sum 

«-l/2 n2 

K{N) := (r32) 

e=o k=o 

by making use of the function /jv(s,t) = N~^C{s,t) exp{NH{s,t)} that appeared in (|4.24p . The corre¬ 
spondence between the two sets of variables {i, k) and (s, t) is given by the simple transformation s = i/N 
and t = k/N. We denote this transformation by Tjv- 

Let us define ^*n and introduce 0 < < min{p, (1 — p)/2}. The precise value 

of 5 will be chosen later on. We partition the domain of the summation appearing in the sum of B^^^{n 2 ) 

Ajv = {(^, k)-. ^ = ni/2, k = l,... 712 }, 

into two subsets 

U5,n = {(^, k)eA^-.\l- Cl < + 1, |A: - Cl < <5A + 1}, CV = Aiv\Civ- 

The set Civ is to be considered as a neighborhood of th® “maximum” point of i?^^^(n 2 ). We 

observe that rjv(C,iv) is contained in the neighborhood of {s*,t*) in Dp, i.e., 

~ ^ Bp-. \s — s*\ < 6 + \t — t*\ < 5 + —}, 

while Tpf{Uf is contained in its complement := Dp\W^j^j_. We rewrite (I4.32p as K(N) = 

Ki{ 6 ,N) + K 2 { 6 , N) where 

K,{5,N):= K2{5,N):= ^ i?5'(n2). (4.33) 

(i,k)eUs,N ikk)&uij, 

We aim to prove that the asymptotic behavior of K{N) is given by Ki{5,N), while K 2 { 5 ,N) gives a 
sub-dominant contribution. We start by proving the latter statement. 
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Bound on K2(5, N). Making use of (14.271) . we upper bound 


K 2 { 6 ,N) e^p{NH{e/N,k/N)} . (4.34) 

Defining 

M{ 6 ) := sup H{s,t) > sup H{s,t) , 

|s-s*|>5 

since the values {£/N,k/N) in (I4.34p belong to kk? j_ we can bound H{£/N,k/N) < M{5). We conclude 

°+~N 

that 

K 2 { 6 ,N) < CN^/^e^p{NM{S)}\UlJ < exp{NM{S)}, (4.35) 

which, together with M{ 6 ) < implies that 

ex.p{-NH{s*,t*)}K 2 { 6 ,N) ^0, as N ^ oo. (4.36) 

Let us remark that, besides the condition 0 < (5 < min{p, (1 —p)/2} (which guarantees that (i, k) and (s, t) 
are contained in the domains of Bf^^(n 2 ) and /jv(s,t)), in the previous argument no further condition has 
been imposed on 5. 


Asymptotics of Ki (<5, N) . Here we consider the sum Ki (<5, N) defined in (I4.33p . Choose e > 0 arbitrary 
and small. By continuity of C{s,t), we can choose <5 > 0 small enough so that, for large N, 

C{s*,t*) — e < C{s,t) < C{s*,t*) + e, for all (s,t)€H( 5 ^j_ 

Then using (|4.24p . we obtain 


^ exp|jVJf(L,L)l(i+„(i,), (4.37) 

ii,k)&Us,N ^ ^ ^ ^ 

and a similar lower bound with C{s*,t*)+£ replaced by C{s*,t*) — e. Recalling that Q(s, t) is the Hessian 
matrix of H{s,t), by Taylor expanding up to second order and using that is the maximum, we can 

write 

H{s, t) — H{s* ,t*) < -X • Q{s* ,t*)x. + c(5||x|p/2, for all (s, t) G 

where x = {s — s*,t — t*), and a similar lower bound with c6 replaced by —c5. 

By multiplying (I4.37P by exp[—A77(s*, t*)] and applying the previous inequality, we obtain 

exp[-AH'(s*,C)]Ai(J,A) 

^ C(s ,^) + e ^ exp jyx^ • Q(s*,C)x + c5A||x||V2| , (4.38) 

{£,k)GUs,N ^ ^ 

(where x is computed with s = £/N and t = k/N) and a similar lower bound with C{s*,t*) + e replaced 
with C{s*,t*) — £ and +c5A||x|p replaced with —c5A||x|p. The last step is bounding the sum 

Ki{ 6 ,N) := Y exp Iyx'^ • (Q(s*,C) ± C(fl)x| , 

{£,k)eU5,N ^ ^ 
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Now we can substitute the finite sum in the previous display with the infinite one, since the difference 
is exponentially small [9]. It is known that (as can be seen by extending [9l (3.9.4)] to two-dimensional 
sums) that 


g-j^Aj/(27V) 

jez2 


2'kN 

det(A)^/2 


(1 + 0 ( 1 )), 


Therefore, 


Ki{5,N) 


2-k 

det(Q(s*, t*) ± C(5I)^/2 


fV(l + o(l)). 


From the previous equation, recalling (I4.38p . we obtain 


exp[-iVi?(s*,r)]iFi((5,iV) < 


2-k{C{s*X) + e) , „nii 

det(Q(s*, t*) — c(fl)^+ ’ 


and a similar lower bound with C{s*,t*) + e replaced with C{s*,t*) —s and Q{s*,t*) — c61 with Q(s*, +) -|- 
c6\. Since e is arbitrary, the previous inequality implies 


lim exp[-iVi/(s*,r)]iF(A) 


27rC(5*,+) 
det(Q(s*, t*))^/2 ’ 


which proves the claim. 

Next, we want to generalize the previous result by computing the asymptotic of (I4.2ip in the case 
m 7 ^ 0. We start by rewriting (14.211) in the following fashion: 


oo ni/2n2—m 

Qn^ci{B)^^ \M^ = m\ = Y ^ , 

1=2 i=0 k=0 


where 




bS’<" 2 - ™) nr.i (=T+"2 - i) n"r„‘(,>2 - k - j) 

B'ui’n) n”r„' («2 - i) n” 1 i^+n^-e-k-j)' 

for m = 0,1,..., n 2 - By defining the function F{s, t; m) on 0° given by 


F{s, t; m) = 


l+p 

2p 


p — t 


i+p 


— s — t 


we obtain that 


G^l'l{m;n2) = F { (l-bo(l)) 


^e,k 


I k 




as N ^ oo. Then the proof is obtained from that for m = 0 by replacing C{s,t) by C{s,t)F{s,t)- D 
Remark 4.1 (Bound on b*). Since {s*,t*) E D°, 


p-t* 


p-t* 


i+p 


- s*) + {p - t*) 


< 1, so that b* < 


l+p 

2p 


This will allow us to use in the following the moment generating function (5jv[(&*)^] defined in (I4.18p . 
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Boundary contribution. Lemma 14.41 proves, for any fixed 0 < m < oo, the asymptotic exponential 
growth of I = m] as iV —?• oo. However in formula (|4.16p we need to sum over a 

range of values of m that increases with the volume N. In order to overcome this problem, in the proof 
of Proposition 14.11 we introduce a cut-off in the sum over m (and then send the cut-off to infinity at the 
end). In doing so we need to exclude the contribution arising from for i close to the boundary ni/2. 
This is achieved in the following Lemma. 

Lemma 4.5 (Boundary contribution). For every e > 0 sufficiently small, as N ^ oo 

712 '^l/2 712—771 

^ ^m[Z^^{B)] Y. E (^2 - = m)= o(e^^(^*’**)). 

m=0 £>(i_£):^ k=0 

Proof. By Lemma [4.31 we know that s* < Defining 

:= sup H{s,t), 

{s,t) 


it follows that H.{s*,t*) < Further, we define 

B>\^^{ni,n2,m) = Bf^j^{n2 - m)QM{MM = 
and using (|4.17p and the following bound 


9-2(m+l) /2(m -F 1)\ 

y m -|- 1 ) 


< 



m). 


we obtain 

D[^^{ni,n2,m + 1) ^ 

712 — Ta — k 

D^^'(ni,n 2 ,m) “ ni/2 - f-|-n 2 - m - A: - 1 

^ _ n2 _^ n2 

~ ni/2 — A-|-n 2 — 1 ~ n 2 — 1 

As a consequence, using (I4.27p . 

- m)QM{M,si = m) < H^^’(n 2 )Qiv(Mjv = 0) 

< aN^/‘^ exp{7L(£/Af, k/N)], 


(4.39) 


/ \ m / \712 

since, ( j < ( n 2 -i ) — ® m < n 2 and some a > e. Therefore, using this inequality together 

with Lemma l3.ll we obtain that 


n2 ni/2 rL2—7n 

m=0 £>(!_£) lU fc=0 

'ni/2 n2 

e>{l-e)^ k=0 


□ 
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Now we are finally ready for the proof of Proposition 14.11 We treat first the case in the presence of an 
external field B and then the case without field. 


Proof of Proposition ^. 1\ (a). We fix G {0,..., n 2 } and e > 0 sufficiently small. Using (|4.16p and Lemma 
we write 

OO 

1=2 


where 




(1 112—m 


^ N,e<(l-e)‘. 


m=0 

n2 


i=0 k=0 

,(//)= ^m[Z^^\B)] Y. Y Bl%\n2-m)Q^{M, = m), 

£=0 k=0 

712 712—771 

^Sxi-e):^ = E E E (^2 - = m). 

m=0 £>(i_£)!^ fc=0 


J^(2) 

N (1 — £} 


(4.40) 

(4.41) 

(4.42) 


We analyze the three pieces separately, showing that only the first of them contributes to the exponential 
growth of Qjv [■^mJ^(.B) ]))[^2 q(-^)'^*] • By Lemma WTU and Lemma 141^ 


=2vrC(s*,r) (det Q(s*,r ))-2 exp{NH{s*,t*)} 


X Y ^m[Z^\B)]Q^iM = m)(l + o(l)). 


(4.43) 


m=0 


The expression in (|4.4ip can be rewritten as 


772 


a \ ''^1 

n 2 —m r(JV) 


= E E, 4 Zff(B)l E E 


B\^^{n 2 - m)Q^{Mj.f = m) 

B‘^’(n2)Qiv(M^ = 0) 




Now, by (|4.39p, 


("2 {Mpf=m) 


is uniformly bounded by 1 — <5 for 5 > 0 sufficiently small. 


Bff.(n2-{m-l))QN(MN=m,-l) 

because i < {1 — £)nij2. Using this bound together with Lemma l3 .1 1 yields 


"'l/2 712 n2 

£=Q k=0 m=^+l 


Thus there exist e(^) (with e(^) —)• 0 as /r —>■ oo) such that, uniformly in N, 

exp{iVi/(s*,C)}(l+ o(l)). 

Finally, from Lemma lL5] it results that = o(exp{A^Lf(s*, U)}). Thus, from (14.431) we can 

identify I = H{s*,t*) and 

OO 

J = 27rC(s*,U)(detQ(s*,r))-5 ^ E^[Z^^(5)](6*)™P(M = m). 

771=0 
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We remark the previous expression is well-defined since, from Lemma l3.ll 

OO 

Y,^m[Z^KB)]{b*rF{M = m) < AE[{b*n 

m=0 

which is finite because b* < (see (I4.18P in Lemma [4.11 and Remark 14.ip . 

Proof of Proposition \J7J\ (b). In this case we work with a vanishing external field. Defining tjv := and 
by Taylor expanding (|4.9p around 0, we have 

o(tiv) := I = <^^(1 + 0 ( 1 )), as N ^ OO, 

where C is a constant (whose value actually depends on /3). Thus, 


<3iv[]jQ(iiv)^' I Mjv = m] 


Qivin (1 + I Mjv = m] 

1=2 

= Qn ^ I Mjv = m 

= Qn I = ml (1 + 0(1)). 


1=2 


By writing 


Qt 




m 


(4.44) 


= eCPj:Z2r‘itN)QNip[^'h 

+ eCPj:Z2r‘itN)QNip\^^)Q^ LcPj:Z2r‘itN)(p\^^-QN{p\^y - 1 I = m , 


and using formula (j4.16p . we can rewrite Qjv 01^2 cz(^iv)’^*] = 5'i(iV) -|- S 2 {N) as the sum of two 

contributions, due to the two terms in (14.441) . Now we analyze S'i(A^) and S2{N) as N ^ 00 . First, we 
remark that the sum in the exponential factor converges in this limit. This can be shown by observing 
that r{B) < 1 (see (j3.111) 1. Therefore, calling r* = r(0) and given any e > 0 such that r* -|- e < 1, thanks 
to the convergence of rftjsf) to r*, we have that for all N sufficiently large, 

00 N N 

1=2 1=2 1=2 


where we used the fact that < 1 for I < N and = 0 for all I > N. Since the geometric sum in 
the r.h.s. of the previous display is convergent, the positive series in the l.h.s. is also convergent to some 
positive value Iq. Thus, by inserting the first term of the r.h.s. of (14.441) in (I4.16P and applying bounded 
convergence, we obtain (I4.12P with 

I{t)=IoCy. 

and 

CX> 

J = ^E™[Z^)(0)]P(M = m). 

m=0 

Further, by Lemma 13.11 and the law of total expectation. 


S2{N) < 


ftN)QN(p\^^)Q^ 
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We use the Cauchy-Schwarz inequality to bound 


Qi 


(4.45) 


^ Qn 

and, by Jensen’s inequality and Holder inequality, 
1 ^ Qn 


Ct^j:Z2r\tN){p\^^-QNip\^^)) _ ^ 


1/2 


< [Q 


Ct^VNj:Z2r‘{tN){p\^^-QN{p\^^)) 


1 

y/N 


< 1 + 0 ( 1 ), 


due to the existence of the finite limit of Qjv 




by [HI Lemma 4.3]. There¬ 


fore, also the term in fl4.45p converges to 0, showing that 5*2(A^) gives a vanishing contribution. This 
completes the proof. □ 

4.2 Annealed thermodynamic limits and SLLN: Proof of Theorems 11.71 and 11.81 

Finally, we prove the existence of the thermodynamic limits. 

Proof of Theorem 11.71 The thermodynamic limit of the annealed pressure is given by 

= lim i;^{f3,B) = lim ^ log (Qjv (■^iv(^,-B))). 

Af->oo N^oo iV 

From (14.101) we can rewrite 

oo 

Q^{Z^il3,B)) = X^{B)Al^^\B)Q^[z^^^jB)llciiB) 

and then 


N, 


1=2 


V’(^,^) = logA+ + 

Using Proposition 14.11 we find 


I — p 


2 logH+ + ^lim^ llog{Q,[zl^;^(H)nQ(^)'^']}- 


1=2 


1 — p 


iA(/3,H) = logA+(/3,H) + ^logH+(/3,H) + F(s*,U). 

To prove the existence of the thermodynamic limit of the magnetization, we use Lemma 12.11 and the 
existence of the pressure in the thermodynamic limit. Then, remembering that is the maximum 

point of the function we compute 

d 


+ s*(/3,H)—log [o(/3,H)r2(/3,H)] 

+ f(/3,H)Aiogr(/?,H) 


(4.46) 


In the limit of small external field B by Taylor expanding (|4.9I) one has a{f5,B) = 0{B‘^). Also, from the 
fixed point equations (|4.30p one can check that s*{f3, B) = 0{B^). As a consequence lim5^o+ B) = 0 
for all /3 > 0, and therefore, by the definition in (II.9h . we conclude that there is no phase transition for 
CM^(1,2). □ 
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Proof of Theorem I j.M ’ Again, the SLLN follows immediately from the existence of the annealed pressure 
in the thermodynamic limit and its differentiability with respect to B. See the proof of Theorem 11.21 and 
0 Section 2.2], □ 
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